######################################################
# 
# --- Command log file ---
# PLAXIS 2D 2025.1
# Tutorial - 13
# Potato field moisture content
# Note: Use the PLAXIS Commands runner to execute this
# 
######################################################
_initializerectangular SoilContour 0 0 15 5
_setproperties "Title" "Potato field moisture content" "Company" "Bentley Systems Inc" "Comments" "" "UnitForce" "kN" "UnitLength" "m" "UnitTime" "day" "UnitTemperature" "K" "UnitEnergy" "kJ" "UnitPower" "kW" "WaterWeight" 10 "ReferenceTemperature" 293.15 "LiquidSpecificHeatCapacity" 4181.3 "LiquidThermalConductivity" 0.0006 "LiquidLatentHeat" 334000 "LiquidThermalExpansion" 0.00021 "LiquidTemperature" 293.15 "IceSpecificHeatCapacity" 2108 "IceThermalConductivity" 0.00222 "IceThermalExpansion" 5E-5 "VapourSpecificHeatCapacity" 1930 "VapourThermalConductivity" 2.5E-5 "VapourSpecificGasConstant" 461.5 "UseTemperatureDepWaterPropsTable" False "ModelType" "Plane strain" "ElementType" "15-Noded"
_borehole 0.75
_borehole 2
_soillayer 0
_set Soillayer_1.Zones[0].Top 3.75
_set Soillayer_1.Zones[0].Bottom 3
_soillayer 0
_set Soillayer_2.Zones[0].Bottom 0
_set Soillayer_1.Zones[1].Top 5
_set Soillayer_1.Zones[1].Bottom 3
_set Borehole_2.Head 4.25
_set Borehole_1.Head 4.25
_soilmat "Identification" "Loam" "SoilModel" "Linear Elastic" "Colour" 15262369 "gammaUnsat" 19 "gammaSat" 19 "ERef" 1000 "nu" 0.3   "GroundwaterClassificationType" "Staring" "SubsoilTopsoil" "Topsoil" "GroundwaterTopSoilClassStaring" "Clayey loam (B9)"  "GwUseDefaults" True
_soilmat "Identification" "Sand" "SoilModel" "Linear Elastic" "Colour" 10676870 "gammaUnsat" 20 "gammaSat" 20 "ERef" 10000 "nu" 0.3 "GroundwaterClassificationType" "Staring"    "SubsoilTopsoil" "Subsoil" "GroundwaterSubSoilClassStaring" "Loamy sand (O2)" "GwUseDefaults" True
_setmaterial Soillayer_1 Loam
_setmaterial Soillayer_2 Sand
_mergeequivalents Geometry
_gotomesh
_set (BoundaryLine_1.CoarsenessFactor BoundaryLine_2.CoarsenessFactor BoundaryLine_4.CoarsenessFactor) 0.5
_mesh 0.06
_gotostages
_set InitialPhase.DeformCalcType "Flow only"
_activate (GWFlowBaseBC_9) InitialPhase
_set GWFlowBaseBC_9.Behaviour InitialPhase "Head"
_set GWFlowBaseBC_9.Href InitialPhase 3
_set GroundwaterFlow.BoundaryXMin InitialPhase "Closed"
_set GroundwaterFlow.BoundaryXMax InitialPhase "Closed"
_phase InitialPhase
_set Model.CurrentPhase Phase_1
_set Phase_1.Identification "Transient phase"
_set Phase_1.PorePresCalcType "Transient groundwater flow"
_set Phase_1.TimeInterval 15
_set Phase_1.MaxStepsStored 250
_dischargefunction
_set DischargeFunction_1.Signal "Table"
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_set DischargeFunction_1.Table[1].Time 1
_set DischargeFunction_1.Table[1].DeltaDischarge 0.01
_set DischargeFunction_1.Table[2].Time 2
_set DischargeFunction_1.Table[2].DeltaDischarge 0.03
_set DischargeFunction_1.Table[3].Time 3
_set DischargeFunction_1.Table[4].Time 4
_set DischargeFunction_1.Table[4].DeltaDischarge -0.02
_set DischargeFunction_1.Table[5].Time 5
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_set DischargeFunction_1.Table[6].Time 6
_set DischargeFunction_1.Table[6].DeltaDischarge 0.01
_set DischargeFunction_1.Table[7].Time 7
_set DischargeFunction_1.Table[7].DeltaDischarge 0.01
_set DischargeFunction_1.Table[8].Time 8
_set DischargeFunction_1.Table[9].Time 9
_set DischargeFunction_1.Table[9].DeltaDischarge -0.02
_set DischargeFunction_1.Table[10].Time 10
_set DischargeFunction_1.Table[10].DeltaDischarge -0.02
_set DischargeFunction_1.Table[11].Time 11
_set DischargeFunction_1.Table[11].DeltaDischarge -0.02
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_add DischargeFunction_1.Table
_set DischargeFunction_1.Table[12].Time 12
_set DischargeFunction_1.Table[12].DeltaDischarge -0.01
_set DischargeFunction_1.Table[13].Time 13
_set DischargeFunction_1.Table[13].DeltaDischarge -0.01
_set DischargeFunction_1.Table[14].Time 14
_set DischargeFunction_1.Table[15].Time 15
_activate Precipitation Phase_1
_set Precipitation.TimeDependency Phase_1 "Time dependent"
_set Precipitation.DischargeFunction Phase_1 DischargeFunction_1
_activate Precipitation Phase_1
_selectmeshpoints
/output addcurvepoint "Node" Soil_1_1 1.75 4.31
/output addcurvepoint "Node" Soil_1_1 1.33 3.80
/output addcurvepoint "Node" Soil_1_1 2.34 3.84
/output addcurvepoint "Node" Soil_1_1 2.48 2.96
/output addcurvepoint "Node" Soil_1_1 5.01 3.52
/output addcurvepoint "Node" Soil_1_1 4.99 2.95
/output addcurvepoint "Node" Soil_1_1 8.63 3.23
/output addcurvepoint "Node" Soil_1_1 8.75 2.96
/output addcurvepoint "Node" Soil_1_1 13.31 3.23
/output addcurvepoint "Node" Soil_1_1 13.40 2.95
/output update